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ABSTRACT 

This paper reports results on the statistical analysis of HI turbulence in the Small Magellanic Cloud 
(SMC). We use 21 cm channel maps, obtained with the Australia Telescope Compact Array and the 
Parkes telescope, and analyze the spectrum of observed intensity fluctuations as a function of the velocity 
slice thickness. We confirm predictions by Lazarian & Pogosyan (2000) on the change of the power law 
index and establish the spectra of 3-D density and velocity. The obtained spectral indices, —3.3 and 
—3.4, are slightly more shallow than the predictions for the Kolmogorov spectrum. This contrasts to 
the predictions for the shock-type spectra that are steeper than the Kolmogorov one. The nature of the 
energy injection in the SMC is unclear as no distinct energy injection scales are observed up to the entire 
scale of the SMC. 

Subject headings: ISM: general - ISM: kinematics and dynamics - ISM: structure - galaxies: individual 
(SMC) - turbulence 
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1. INTRODUCTION 

Many observations in the past decade have challenged 
the traditional picture of the interstellar medium (ISM). 
Instead of a two-level hierarchical system, consisting of 
clouds uniformly dispersed in the intercloud medium (the 
so called 'standard cloud' model, Spitzer 1978), the ISM 
shows an astonishing inhomogeneity, with many levels of 
hierarchy. In order to consider real density functions in 
physical processes, a better understanding of the inventory 
and topology of the ISM is essential, as well as of the pro- 
cesses responsible for their creation. Having an extremely 
gas-rich ISM, dwarf irregular galaxies are particularly suit- 
able for such studies. Here we discuss the inventory of the 
ISM in the Small Magellanic Cloud (SMC), using the spa- 
tial power spectrum, and point to several processes that 
may be involved in the sculpturing of its ISM. 

The SMC is a nearby 1 , extremely gas-rich, dwarf irreg- 
ular galaxy. Taking a part in an interacting system of 
galaxies (with other members being the Large Magellanic 
Cloud, LMC, and our Galaxy), the SMC's morphology, 
dynamics and evolution is very complex. Being one of our 
closest neighbors, the SMC is attractive for various astro- 
physical aspects. The inventory of the cool ISM in the 
SMC was recently studied using the high resolution radio 
observations of neutral hydrogen (HI) (Stanimirovic et al. 
1999). Also, the dust properties of the SMC were investi- 
gated in detail, and the relationship between the cool gas 
and dust, using both HI and infrared (IR) observations 
(Stanimirovic et al. 2000). 

As a dwarf, irregular galaxy, the SMC is different from 
our own Galaxy in many respects: its interstellar radiation 
field (ISRF) is more than four times stronger (Lequeux 
1989), its heavy element abundance is almost ten times 
lower (Sauvage & Vigroux 1991), its dust grains are on 
average smaller by « 30% (Rodrigues et al. 1997) and sig- 



nificantly hotter (Stanimirovic et al. 2000), and the cool 
atomic phase of HI is only half as abundant (Dickey et al. 
2000). In spite of all the above, the SMC appears to have 
similar interstellar turbulence properties, in having similar 
power law index as the Galaxy for the 2-D spatial power 
spectrum of its HI and dust column density distributions 
(Stanimirovic et al. 2000). 

What do the 2-D spatial power spectra of the intensity 
fluctuations mean? Relating the fluctuations of intensity 
in PPV (position-position- velocity) data cubes and the un- 
derlying 3-D velocity and density statistics is a problem 
that has been recently addressed in Lazarian & Pogosyan 
(2000), where it was shown that changing the thickness 
of the velocity slice it is possible to recover both spectra 
of turbulent velocity and density. The SMC provides an 
ideal testing ground for this theory and we apply theo- 
retical predictions to the data. Elsewhere we plan to ap- 
ply alternative tools for turbulence studies, e.g. Principal 
Component Analysis (PCA) (Heyer & Shloerb 1997; Brunt 
& Heyer 2001), the A-variance (Stutzki et al. 1998) and 
the spectral correlation function (Rosolowsky et al. 1999) 
(see review by Lazarian 1999) to the SMC data. 

The structure of this paper is organized as follows. In 
Section 2 we describe briefly the HI observations of the 
SMC. In Section 3 we review previous results of the sta- 
tistical investigation of HI in the SMC, using the spatial 
power spectrum, and investigate the influence of velocity 
fluctuations on the intensity statistics. Comparison with 
earlier work and discussion on the origin of the turbulence 
in the SMC, are given in Section 4. 

2. HI DATA 

The small-scale HI structure in the SMC was observed 
with the Australia Telescope Compact Array (ATCA), a 
radio interferometer, in a mosaicing mode (Staveley-Smith 
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et al. 1997). Observations of the same area were obtained 
also with the 64-m Parkes telescope, in order to map the 
distribution of large-scale features. Both sets of obser- 
vations were then combined (Stanimirovic et al. 1999), 
resulting in the final HI data cube with the angular reso- 
lution of 98 arcsec, the velocity resolution of 1.65 km s _1 , 
and the 1-er brightness-temperature sensitivity of 1.3 K, to 
the continuous range of spatial scales between 30 pc and 
4 kpc. The velocity range covered with these observations 
is 90 to 215 km s _1 . For details about the ATCA and 
Parkes observations, data processing and data combina- 
tion (short-spacings correction) see Staveley-Smith et al. 
(1997) and Stanimirovic et al. (1999). 

3. STATISTICAL INVESTIGATION 

The 2-D spatial power spectrum, P(k), of the HI emis- 
sion fluctuations (I) in the SMC was first derived in Sta- 
nimirovic et al. (1999). This was the first such study for 
an entire galaxy. The power spectrum is defined as the 
Fourier transform of the autocorrelation function of the 
HI emission fluctuations: 

P(k)= J y(/(x)/(x'))e- tLk (iL,L = x-x' (1) 

with k being the spatial frequency, measured in wave- 
lengths, and L being the distance between two points. 
To derive P(k) the channel maps (approximately 5 
km s" 1 wide) were Fourier transformed. The average 
value of the square of the modulus of the transform, 
(jft 2 + 3 2 ), was then measured in 18 annuli of equal width 
in log y/u 2 + v 2 (u and v are the coordinates in the Fourier 
plane, being measured in wavelengths) assuming radial 
symmetry. It was shown that the 2-D power spectra were 
remarkably well fitted by a power law, P(k) oc k 1 , over the 
continuous range of spatial scales <~ 30 pc - 4 kpc, and over 
the velocity range 110 - 200 km s _1 (several noisy chan- 
nels, without significant line emission, were excluded from 
the analyses). The power- law index was remarkably con- 
stant over the given velocity range, and the average slope 
was estimated with (7) = —3.04 ± 0.02. No change of the 
power law slope was seen on either large or small-scale end 
of the spectrum. 

However, an examination of the power spectrum of the 
HI column density distribution, obtained after integrating 
along the whole velocity range, showed a change of the 
power law slope, with 7 being equal —3.31 ± 0.01 (Stan- 
imirovic et al. 2000). 

The power-law behavior of the spatial power spectrum 
was seen before in the case of our Galaxy, for several ISM 
tracers (Crovisier & Dickey 1983; Green 1993; Gautier 
et al. 1992; Waller & Varosi 1997; Schlegel et al. 1998; 
Deshpande et al. 2000), and very recently in the case of 
the LMC (Elmegreen et al. 2000). The integrated HI col- 
umn density distribution of the SMC, has a steeper spatial 
power spectrum (7 = — 3.3±0.01), which is consistent with 
the spectrum of the dust column density (7 = —3.1 ± 0.2) 
(Stanimirovic et al. 2000). In several cases, a steepening 
of the slope for higher spatial frequencies was found (Cro- 
visier & Dickey 1983; Waller & Varosi 1997; Elmegreen 
et al. 2000). A change of power-law slope when examining 
channel maps and the integrated intensity along the whole 
velocity range, noticed in Stanimirovic et al. (2000), was 
not found previously. 



3.1. Velocity modification of the HI power spectrum 

The 2-D spatial power spectrum of the intensity (or 
emissivity) fluctuations provides important information 
about the underlying 3-D statistical properties of the ISM. 
However, as the spectral-line data cubes have two spatial, 
and one velocity axis, both 3-D density and velocity fluc- 
tuations can contribute to the 2-D statistics. Indeed, due 
to the velocity fluctuations, two clumps along the same 
line of sight at different distances may appear in the same 
velocity channel, hence doubling the measured intensity. 
In order to relate the 2-D statistics to the underlying 3-D 
statistics of the ISM, one must disentangle density from 
velocity influence to the power spectrum. 

The velocity statistics is probably the most interesting 
and significant, partly because it directly represents the 
dynamics of the media, and partly because other methods 
of turbulence study, i.e. scintillation technique (Spanglcr 
1999), fail to deliver this part of information. The prob- 
lem of disentangling velocity and density fluctuations was 
addressed in Lazarian & Pogosyan (2000) (see also Lazar- 
ian 1999 for a review of the subject). They found that 
the relative influence of velocity and density fluctuations 
on the intensity fluctuations depends on the spectral in- 
dex of density fluctuations and on the thickness of channel 
maps. The latter is easy to understand qualitatively, as it 
is clear that line integration should decrease the influence 
of velocity fluctuations. The former is less intuitive but 
Lazarian & Pogosyan (2000) showed that if the 3-D den- 
sity spectrum is P s oc k n , two distinct regimes are present 
when: (a) n > —3, and (b) n < —3. 

Analytical results of Lazarian & Pogosyan (2000) study 
are summarized in Table 1. It is easy to see that in both 
cases (a) and (b) the power law index gradually steepens 
with the increase of velocity slice thickness. In the thickest 
velocity slices the velocity information is averaged out and 
it is natural that we get the density spectral index n. The 
velocity is the most important in thin slices, and if the 3-D 
velocity power spectrum is k~ 3 ~ m , then the index m can 
be found from thin slices. Note that the notion of thin and 
thick slices depends on a turbulence scale under study and 
the slice can be thick for small scale turbulent fluctuations 
and thin for large scale ones. The formal criterion for the 
slice to be thin is that the dispersion of turbulent velocities 
on the scale studied should be less than the velocity slice 
thickness. The slice is thick otherwise. When velocity in- 
formation is averaged out (for a more precise definition see 
Lazarian & Pogosyan 2000) the slice gets very thick and 
the velocity fluctuations get independent of density. 

Predictions of Lazarian & Pogosyan (2000) were suc- 
cessfully tested using numerical MHD data in Vazquez- 
Scmadeni et al. (2000) and thus we can confidently apply 
the theory to observations. 

To test the predictions by Lazarian & Pogosyan (2000) 
in the case of the SMC, we have determined the 2-D power 
spectrum slope, (7) , while varying the velocity slice thick- 
ness from <~ 2 km s _1 to ~ 100 km s _1 . A significant 
variation of (7) was found, shown in Fig. 1, consistent 
with the predictions — (7) gradually decreases with an 
increase of velocity slice thickness. The thickest velocity 
slice gives n s=a —3.3, which corresponds to the 3-D density 
power index (see Table 1). Using thin slices, however, we 
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can find the slope of velocity fluctuations to be m w 0.4, 
which means that the 3-D velocity spectrum index is —3.4, 
which is very close to the 3-D density spectral index. The 
transition point between density and velocity dominated 
regimes is equal to the velocity dispersion on the scale of 
the whole SMC (~ 4 kpc), which is ~ 22 km s" 1 (Stan- 
imirovic et al., in preparation). 

4. DISCUSSION 

As the spatial power spectrum shows the importance of 
structure on various spatial scales, its power-law behav- 
ior suggests the hierarchical structure organization in the 
ISM, without preferred spatial scales. This phenomenon is 
usually ascribed to the interstellar turbulence (Scalo 1987; 
Elmcgrecn 2000). However, without velocity information 
one can always wonder whether we deal with a static struc- 
ture or a real turbulence. Indeed, a distribution of sizes 
of sand grains on a beach also follows a power-law, but no 
one would call this "turbulence" . The velocity information 
changes the picture dramatically. Hence, the extreme im- 
portance of the techniques which relate the observed 2-D 
power spectrum with the underlying 3-D statistics of both 
density and velocity. Here, we have tested the theoreti- 
cal predictions for such a technique (Lazarian & Pogosyan 
2000), and as a result proved, for the first time, the pres- 
ence of an active turbulence in the SMC. 

In view of theoretical results in Lazarian & Pogosyan 
(2000) it is now appropriate to reanalyze all the earlier 
data. These data were obtained without much concern 
about the thickness of velocity slices. Therefore the ob- 
served variations of the power index can be due to transi- 
tions from "thin" to "thick" and to "very thick" slices. In 
the case of Green's (1993) data additional complications 
are related to a divergent line of sight geometry. A more 
detailed discussion of the available data will be given else- 
where. We note that although the description in terms of 
power spectra is common in hydrodynamics and the MHD 
theory, it has certain limitations, as discussed in Lazarian 
(1999). For example, the power spectrum analysis does not 
include information about the phase distribution, dealing 
only with the modulus of the Fourier transform, nor it con- 
tains information about the structure connectivity (Scalo 
1987). Other methods hence, should be used as comple- 
mentary statistical descriptors. Power spectrum as it is 
can provide us with an important insight of what kind of 
turbulence we deal with, e.g. distinguish the turbulence 
originating from shock waves from the hydrodynamic tur- 
bulence. 

The attempts to test Lazarian & Pogosyan (2000) theory 
were made recently in Elmcgrcen et al. (2000) using the 
HI observations of the LMC. In agreement with theoreti- 
cal predictions, the steepening of spectrum was observed 
for high spatial frequencies. The puzzling thing discovered 
by Elmegreen et al. (2000) was the flattening of the spec- 
tra for velocity-integrated intensity, which was interpreted 
as an effect of the finite LMC disk thickness. This is an 
interesting explanation which entails that the LMC spec- 
trum at the scales larger than 100 pc becomes essentially 
two dimensional. Our study has not noticed a systematic 
change of the velocity integrated power spectrum at large 

2 This would correspond to the spectrum of shock waves. To avoid possible confusion we point out that we talk about power spectrum which 
differs from the energy spectrum by k 2 . The Kolmogorov turbulence corresponds to the power spectrum of —11/3. 



spatial scales. This may reflect the fact that the SMC, 
unlike the LMC, is essentially a 3-D entity. 

Another approach in relating the 2-D with the 3-D 
statistics in the case of the SMC was presented in Goldman 
(2000). There it is assumed that the density fluctuations 
are a "passive scalar" , being driven by the velocity fluctu- 
ations, and hence having the same power spectrum. If we 
accept that the intensity fluctuations are due to the den- 
sity fluctuations, then the corresponding spectral index q 
of intensity fluctuations in a 2-D slice can then be related 
to the 3-D density spectrum index n as 

q = n+l , (2) 

which for the SMC data produces n w — 4 2 . However, we 
note that the data used in Goldman (2000) are not in the 
real space (xyz) for which his treatment would be correct, 
but are in the velocity space (xyv). In this situation the 
Lazarian & Pogosyan (2000) treatment is appropriate and 
it provides a different result, namely, velocity index w — 3.4 
and density index w —3.3. Eq. (2) also predicts that the 
difference in the power slope between thin and thick slices 
is equal to 1. This is inconsistent with Fig. 1. We also 
note that for the Kolmogorov spectrum the predictions in 
Lazarian & Pogosyan (2000) (see Table 1) coincide with n 
calculated using Eq. (2), but this correspondence is acci- 
dental. 

An interesting application of the power spectrum of 
HI opacity fluctuations was made by Dcshpandc (2000), 
in order to explain the long-standing puzzle of the tiny- 
scale structure in HI (Heiles 1997). Assuming a single 
power spectrum of the opacity fluctuations, with a slope 
of 2.75 over the range of ~ 0.02 pc to <~ 4 pc, Deshpande 
(2000) obtained opacities consistent with the observations 
of small-scale HI structure in Deshpande et al. (2000). 
This is very encouraging and requires re-interpretation of 
previous observations of the small-scale structure in a sim- 
ilar way. A preliminary investigation of the small-scale 
structure found so far by Heiles (2000) suggested though a 
more complex structure function, with a significant change 
of slope for scales smaller than 0.01 pc. 

What can drive the turbulence in the SMC? The natural 
assumption would be that it is due to the stirring of the 
ISM produced by a large number of expanding shells found 
in the SMC. The shell sizes range from ~ 30 pc to ~ 2 kpc. 
Hence, one scenario could be that the largest shells drive 
the turbulent cascade down to smallest observed scales. 
However, processes like shell fragmentation and/or shell 
propagation from the smaller scales (bottom-up scheme, 
see Scalo 1987), may play significant role too. The main 
problem in pinning down the exact mechanisms is that, so 
far, we have not observed changes in the power spectrum, 
at any scale up to the entire size of the SMC, which would 
be indicative of energy injection. An alternative explana- 
tion was suggested in Goldman (2000), whereby the large 
scale turbulence is induced by instabilities in the large- 
scale flows during the last SMC-LMC encounter. Future 
comparison with simulations of different types of turbulent 
cascades are essential to resolve this question. 

5. CONCLUSIONS 
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We have successfully tested predictions of the Lazarian 
& Pogosyan (2000) study on the change of slope of the in- 
tensity fluctuation spectrum with the velocity slice thick- 
ness. The SMC spectrum appears due to active turbulent 



motions rather than just a static hierarchical structure. 
We found the index for the 3-D velocity spectrum to be 
—3.4, while for the density spectrum to be —3.3. 
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Table 1 

a summary of analytical results for the 2-d power spectrum of intensity fluctuations, p(k), derived in 
Lazarian & Pogosyan (2000) for two distinct regimes of the 3-D density fluctuations. The index of the 3-D 

POWER SPECTRUM OF VELOCITY FLUCTUATIONS IS —3 — m. 
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Fig. 1. — The variation of the 2-D HI power spectrum slope (7} with the velocity slice thickness Av. The dot-dashed lines distinguish thin, 
thick and very thick slice regimes. 



